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ABSTRACT 

The Crab pulsar is a quite young famous pulsar which radiates multi-wavelength pulsed photons. The latest 
detection of GeV and TeV pulsed emission with unprecedented signal-to-noise ratio, supplied by the powerful 
telescopes: Fermi, MAGIC and VERITAS, challenges the current popular pulsar models, which can be a 
valuable discriminator to justify the pulsar high-energy-emission models. 

Our work is divided into two steps. First of all, taking reasonable parameters (the magnetic inclination angle 
a = 45° and the view angle Q — 63°), we use the latest high-energy data to calculate radio. X-ray, 7-ray and 
TeV light curves from a geometric view to obtain some crucial information on emission locations. Secondly, 
we calculate the phase-averaged spectrum and phase-resolved spectra for the Crab pulsar and take a theoretical 
justification from a physical view for the emission properies as found in the first step. It is found that a Gaussian 
emissivity distribution with the peak emission near the null charge surface in the so-called annular gap region 
gives the best modeled light curves. The pulsed emission of radio. X-ray, 7-ray and TeV are mainly generated 
from the emission of primary particles or secondary particles with different emission mechanisms in the nearly 
similar region of the annular gap located in the only one magnetic pole, which leads to the nearly "phase- 
aligned" multi-wavelength light curves. The emission of peak 1 (PI) and peak 2 (P2) is originated from the 
annular gap region near the null charge surface, while the emission of bridge is mainly originated from the core 
gap region. 

The charged particles cannot corotate with the pulsar and escape from the magnetosphere, which determines 
the original flowing primary particles. The acceleration electric field and potential in the annular gap and core 
gap are huge enough in the several tens of neutron star radii. Thus the primary particles are accelerated to ultra- 
relativistic energies, and produce numerous secondary particles (pairs) in the inner region of the annular gap 
and core gap. We emphasize that there are mainly two types of pairs, i.e., one is curvature-radiation induced 
(CR-induced), and the other is inverse-Compton-scattering induced (ICS-induced). The phase-averaged spec- 
trum and phase-resolved spectra from soft X-ray to TeV band are produced by four components: synchrotron 
radiation from CR-induced and ICS-induced pairs dominates the X-ray band to soft 7-ray band (lOOeV to 
lOMeV); curvature radiation and synchrotron radiation from the primary particles mainly contribute to 7-ray 
band (lOMeV to ^ 20GeV); ICS from the pairs significantly contributes to the TeV 7-ray band (^ 20 GeV to 
400 GeV). 

The multi-wavelength pulsed emission from the Crab pulsar can be well modeled with the annular gap and 
core gap model. To distinguish our single magnetic pole model from two-pole models, the convincing values 
of the magnetic inclination angle and the viewing angle will play a key role. 

Subject headings: gammaiays: stars -pulsars: general -pulsars: individual (Crab: PSRB053lH-21)-radiation 
mechanisms: non-thermal - X-rays: individual (PSR B0531H-21) - acceleration of particles 



L INTRODUCTION 

The Crab pulsar (PSR B0531H-21 or PSR J0534H-2200) is 
the second most energetic pulsar to date with a spin-down lu- 
minosity E'lot — 4.6 X lO'^^crgs"^. It is at a distance of 
d ^ 2kpc, and located in the Crab Nebulae which is a center- 
filled remnant of a supernova discovered by the Chinese "as- 
tronomers" in 1054 AD. The Crab pulsar, a very young pul- 
sar with a characteristic age r = 1240 yr and a characteristic 
magnetic fie ld of Bq — 3.715 x 10 ^^ G, has a spin period of 
P = 33 ms ([Manchester et al.ll2005b . 

The Crab pulsar radiates multi-wavelength pulsed emission 
fromradio (10^^ eV) to 7-ray (up to TeV) band, especially the 
new 7-ray results of the sensitive Large Area Telescope (LAT) 
on board the Fermi Gamma-ray Space Telescope (Fermi) and 
a discovery of TeV pulsed emission by the powerful VERI- 
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TAS array of atmospheric Cherenkov telescopes. Fermi LAT 
presented the high-quality 7-ray (lOOMeV to 20 GeV) light 
curves and spectral data using 8 months of survey data, and 
predicted an exponen tial power-law sp ectra with cut-off en- 
ergies of a few GeV (lAbdo et al.ll20Toh . To verify the expo- 
nential power-law cut-off spectra, the 25 — 100 GeV pulsed 
emission from the Crab pulsarhas been precisely measured by 
MAGIC telescope (lAleksic et al.ll2011l) . It is shown that the 
observed cut-off spectr um has a large deviation from the in- 
ferred exponential one. [VERITAS CoUaboration etall (|20TT1) 
reported that the pulsed emission above 1 00 GeV from the 
Crab pulsar has been detected by the VERITAS, and showed 
their light curves and spectral data with excellent signal-to- 
noise ratio. They declared that current popular pulsar models 
(e.g. outer gap and slot gap models) cannot explain the de- 
tection, and the observation might not be explained by the 
curvature radiation as the origin of the observed emission 
above 100 GeV. These findings enable us to obtain consid- 
erable insights of the magnetosphere physics, e.g., accelera- 
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tion electric field, emission region and the relevant emission 
mechanisms. The multi-wavelength phase-averaged spectrum 
and phase-resolved spectra can discriminate the various pul- 
sar non-thermal emission models. 

There are four major physical or geometrical magneto- 
spheric models which have previously been proposed to 
explain pulsed 7-ray emission of pulsars: the polar cap 
model ('Daughertv & Harding 1 11994 lig QeV the outer gap 
model (Cheng etal. 1986a b: Romani & Yadigaroglu 1995; 



Zhang & Cheng 199 71: ICheng et al. 2000; Zhang et al. 2004, 
2007i:",Hirotani " 2008riTang et al. 2008: ,Lin & Zhang 2009) . 



the two-pole caustic ( TPC) model or the slot gap model 
(iDvks & Ru dald '20Q1 iMuslimov & Harding! |2^ llOOl 
Harding et a l. 2008) , and the annula r gap model ( lOiao et alj 
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2004allbll2007l:lDu et al.ll20Toll201 lb . The distinguishing fea- 



tures of these pulsar models are different acceleration electric 
field region for primary particles and rel evant emission mech- 
anisms to radiate high energy photons jDu et al.ll201 II) . One 
of the key discrepancy of the mentioned emission models is 
the two important geometry parameters: the magnetic incli- 
n ation angle a and th e view angle C^. 

ICheng et al.l (120001) used a 3D single-pole outer gap model 
to present incipient results of light curves, phase-averaged 
spectrum and ph ase-resolved spec tra for the Crab pulsar 
Eight years later, iTang et aTl (l2008h improved the outer gap 
model and used a modified outer gap model considering the 
emission from both poles to calculate light curves, phase- 
averaged spectrum and phase-resolved spectra from lOOeV 
to 10 GeV. In the case of the larger viewing angle C = 78° 
or C = 83 ° and intermediate inclination angle a = 45°, 
IZhang&Lil(l2009h :lLi & Zhang (2010) improved the 3D two- 
pole outer gap model dTang et al.l l2008) to present their best 
results of light curves, phase-averaged spectrum and GeV 
phase-resolved spectra. They also used physical emissivi- 
ties to calculate the light curves, showed the best-fit phase- 
averaged spectrum using emission components from the both 
magnetic poles, and indicated that I CS from pairs mainly con- 
ti-ibuted to the 10"* - 10^ ke V band. iHarding et all (l2008h used 
a 3D slot gap model developed from the TPC model to cal- 
culate optical-to-7-ray light curves, phase-averaged spectrum 
and phase-resolved spectra for the Crab puls ar with as s uming 
a broken power-law pair energy spectrum. iHirotanil (l2008h 
demonstrated that the slot gap model reproduces at most 20% 
of the observed GeV fluxes owing to the small trans-field 
thickness. 

To well study pulsed emission from pulsars, there may be 
two methods: (1) one can use a physical model with a realistic 
accelerating electric field to directly calculate the light curves 
and spectra; (2) one can also utilize a reasonable assumption 
of numerical emissivity to calculate light curves from a geo- 
metric point of view; then obtain some valuable information 
of radiation locations, finally calculate the spectra and take 
a consistent theoretical justification from a physical point of 
view. In conventional cases, some scientists chose the first 
method. Here we choose the second way to study pulsed 
emission from the Crab pulsar. In this paper, we study multi- 
wavelength light curves, phase-averaged spectrum and phase- 
resolved spectra of the Crab pulsar in the annular gap model. 
In § 2, we briefly introduce the annular gap and core gap, pair 
production and calculate the acceleration potential in the an- 
nular gap. In § 3, we model the multi-wavelength light curves 
using the annular gap model together with a core gap. The 
radio emission region is identified and the reason for phase- 
aligned peaks of multi-wavelength light curves is explained. 
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Fig. 1. — The shapes of the annular gap (AG) and core gap (CG) in the 
polai' cap region for the crab pulsar with an incHnation angle of a = 45°. 
The core gap region is defined by the footpoints of critical field lines on the 
polar cap, and the annular gap region is between the footpoints of critical field 
lines and last open field lines. The annular region becomes more symmetric 
if OL is smaller 

To model spectra of the Crab pulsar, we also calculate the 
multi-wavelength phase-averaged and phase-resolved spectra 
of radiation from both primary particles and pairs. Finally, 
conclusions and discussions are presented in § 4. 

2. THE ANNULAR GAP AND CORE GAP 

2.1. Formation of the Annular Gap and the Core Gap 

As noted bv lDu et aTl (1201 lb . the open-field-line region of a 
pulsar magnetosphere is divided into two isolated parts by the 
critical field lines which denote a set of special field lines that 
satisfy the condition of O • B = at the light cylinder The 
core region around the magnetic axis is defined by the criti- 
cal field lines, and the annular region is located between the 
critical field lines and the last open field lines (see Figure[T]i. 
The width of the annular polar region is anti-correlated with 
the pulsar period, it is therefore larger for pulsars with smaller 
spin periods. The annular acceleration potential is negligible 
for older long-period pulsars, but very important for pulsars 
with a small period, e.g., millisecond pulsars and young pul- 
sars. The acceleration electric field extends from the pulsar 
surface to the null charge surface or even beyond it. The an- 
nular gap has a sufficient thickness of trans-field lines and a 
wide altitude range for particle acceleration. In the annular 
gap model, the high energy em ission is g enerated in the vicin- 
ity of the null charge surfac e dDu et alj 
a wide 7-ray emission beam ( lOiao et alj \ 



201c ). This leads to 



20071) . The radiation 



components from both the core gap and the annular gap can be 
observed simultaneously by one observer ( lOiao et al.ll2004bh 
if the inclination angle and the viewing angle are suitable. 

2.2. Acceleration Potential in the Annular Gap 

Here we will explore the formation mechanism of the ac- 
celeration electric field in the annular gap. To give a simpli- 
fied picture for the powerful acceleration electric field along 
an open field line, we attempt to derive ID continuous solu- 
tion for the acceleration potential. Then we use this realistic 
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Fig. 2. — Acceleration potential on an open field line with a magnetic az- 
imuthal \I' = 0° in the annulai' gap. The locations of maximum potential and 
null charge surface are marked. The maximum potential drop on this field 
fine is comparable to the one generated by the unipolar effect. 

acceleration field to obtain the Lorentz factor 7p of the pri- 
mary particle which is from the balance of the acceleration 
and curvature radiation reaction (see Equation (|5]l in §2.3). 
The 7p is not a true maximum Lorentz factor but a valuable 
criterion, it can be regarded as an upper limit of the real max- 
imum Lorentz factor of the primary particles which shape the 
obsered spectra. 

Under the assumption of a fully char ge-separated magne - 
tosphere, the outer gap can be formed (ICheng et alJll986"ah . 
However, if abundant pairs are produced, this lead to a large 
neutral charge component in the magnetosphere. In this case, 
the acceleration mechanism of particles becomes different 
from the case of outer gap model. 

The pulsar magnetosphere is filled with charge-separated 
pair plasma, and the charged particles can not co-rotate with 
the neutron star near the light cylinder and must escape from 
the magnetosphere. This is the generation mechanism for ac- 
celeration electric field. The annular gap and the core gap 
simultaneously export charged particles with opposite sign, 
which can lead to the circuit closure in the whole magneto- 
sphere. The parallel acceleration electric field (E0 in the 
annular gap and core gap regions are opposite. As a result, 
vanishes at the boundary (the critical field lines) between 
the annular and the core regions and also vanishes along the 
closed field lines. The positive and the negative charges are 
accelerated from the core and the annular regions, respec- 
tively. 

To unveil the acceleration potential in the annular gap re- 
gion, we now consider a tiny magnetic tube embedded in an 
open field line. We assume that the particles flow out of the 
co-rotating magnetosphere at a radial distance about rout ~ 
i?LC = 1-57 X 10^ km, and that the charge density of flowing- 
out particles pb(?'out) is equa l to the local Goldreich-J ulian 
(GJ) charge density pgj (rout) (iGoldreich & Julianlll969h . For 
any altitudes r < rout, Ph{r) < Pgjir)- The acceleration 
electric field therefore exists along the field line, and cannot 
vanish until approaching the altitude of rout- 



In a static dipole magnetic field configuration, the field 
components can be described as = 2^cos6> ^^ -g^ _ 

^^'f ^ ng, here 9 is the zenith angle in magnetic coordinate, 
and Bq is the surface magnetic field. Thus the magnetic field 
strength at a altitude r is B{r) = MoRi VScos^ e+i^ 

In the co-rotating frame, the equation for acceleration po- 
tential 5" is 



-47r(pb-Pgj). 



(1) 



Although the same Possion equation for acceleration field 
is used by all of pulsar emission models (some models also 
considered the relativistic effect), the key discrepancy among 
these models is how to obtain and explain the difference 
between flowing charge density and the local GJ density 
(ph — Pgi)- Normally it is assumed that pb — Pgj ~ at 
the star surface. This leads to the quite different results of the 
acceleration electric field. 

Using the conservation laws of the particle number and 
magnetic flux in the magnetic flux tube, the difference be- 
tween the flowing charge density and local GJ charge density 
at the altitude r can be written as 

ilB(r) 

Ph (r) - Pgj{r) = (cos (out - COS C) , (2) 

where ft = 27r/P is the angular velocity, P is the rotation 
period, and C and (^out are the angle between the rotational 
axis and the B field direction at r and rout, respectively. Wang 
et al. (2006) found 



cos C = cos a cos 8^ — sin a sin 8^ cos "0, 



(3) 



where t/j and 8^ are the azimuthal angle and the tangent angle 
(half beam angle) in the magnetic field coordinate, respec- 
tively. Combine Equations ([ij, 0, and (|3]l, we then obtain 
one-dimensional two-order differential equation for the accel- 
eration potential, i.e. 



s' d9 



nBoR^ 



VScos^ 8 + l(cosCout— cosC)s9 , 

(4) 




d9^ Sg atl cr-^ 
where Sg, Sg are given by 



" dSg 

r = Rc sin^ 9, 

here i?o is a field line constant, which denotes for the maxi- 
mum length of a certain point on a field line. Then substituting 
= 2-ta'n^^9 jOiao & LinllI998l) into Equations © and 
(lU, we can solve the Equation ^ , and achieve the 1-D so- 
lution for the electric potential 5* along a magnetic filed line 
with a magnetic azimuthal of ■</; = 0° for the Crab pulsar, as 
shown in Figure |2l The maximum potential drop on this field 
line is comparable to the potential 

Ay = = 1.0 X io^^B,2RlP-\ 

which is generated by the unipolar effect. The acceleration 
potential is quite huge in the inner region of annular gap, 
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and the primary particles are therefore accelerated to ultra- 
relativistic energy with large Lorentz factors of 7 ^ 10^—10^. 
Simultaneously, the accelerated primary particles emit abun- 
dant 7-ray photons through ICS and CR process, then dense 



pairs are generated via 7— B (photon magnetic absorption) 
process. 

2.3. Pair Production 

Since ICS and curvature radiation are two effective radia- 
tion mechanisms to generate high energy 7-ray photons, 
pairs can be generated by these ICS and CR photons emit- 
ted from the accelerated primary particles in both annular gap 
and core gap. Thus three gap modes exist for pair produc- 
tion, namely, CR ga p, thermal ICS gap and resonant ICS gap 
(IZhang et al.lll997allbh . which will be briefly introduced be- 
low. 

In t he traditional inner gap model (iRuderman & Sutherland! 
11975b . 7 — _B process plays a very important role, two con- 
ditions should be satisfied at the same time for pair produc- 
tion: (1) to produce enough high energy 7-ray photons, a 
strong enough potential drop should be reached; (2) for pair 
production, the energy component of 7-ray photons perpen- 
dicular to the magnetic field must satisfy the condition of 

The accelerated particles are assumed to flow along a field 
line in a quasi-steady state. Using the calculated acceleration 
electric field, we can obtain the Lorentz factor 7p of the pri- 
mary particle from the curvature radiation reaction 



7p = ( 



2e 



= 2.36 X 10 V -E^ 



(5) 



where e is the charge of an electron, pr is the curvature radius 
in units of 10^ cm and g is the acceleration electric field in 
units of lO'^Vcm^^ 

In 7— B process, the conditions for pair production are that 
the mean free path of 7-ray photon in strong magnetic field is 
equal to the gap he ight, / h. The mean free path of 7-ray 
photon is given by (lErberill966h 



4.4 h Be ^ 4 
e'^ /hcrrieC B± 3x ' 



(6) 



where = 4.414 x 10^'^ G is the critical magnetic field, H is 
the reduced Planck's constant. 
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X 
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znipc' Br 



E-y B± 

2mpC^ Be 



(7) 



and B± is the magnetic field perpendicular to the moving di- 
rection of 7 photons, which can be expressed as (RS75) 



-B : 



I 



-B. 



(8) 



Here Z sa ft, is the condition for gap sparks (pair production) to 
take place, p is curvature radius of a spot on a magnetic field 
line. For a dipole magnetic configuration, it can be estimated 
as 



-iXRc/ny/' 



(9) 



( IZhang et al.lll997bl) if the spot position is near the neutron 
star surface. Here A is a parameter to show the field lines, A = 
1 corresponding to the last open field fine. The characteristic 



energy from the curvature radiation process can be written as 



E^ ^cr — ^~ o 
Zp 



Then the CR gap height /icr is 



/iCR^5xl03p3/7s-V7„2/7 



(10) 



(11) 



(IZhang etal.lll997bb . 

The so-called CR-mode gap implicates that the pair pro- 
duction cascades are dominated by the 7-ray photons em- 
anated from curvature radiation process of the accelerated 
primary particles, this gap is somehow like the RS gap 
(Ruderman & Sutherland 1975). The resonant ICS gap is 
formed from the resonant scattering of soft photons, which is 
a quantum effect with a large electron scattering cross section 
estimated as cr '--^ IO^ctt (c tt is the Thompson cross section) 
in strong magnetic fields. IZhang et"al] (Il997bl) obtained the 
gap height (/iics) of the resonant ICS mode 



1.1 X 10^P^/^Sf2V6^^ cm 



(12) 



and the Lorentz factor (72,10s) of pairs for the resonant ICS 
mode 



72, 



890P-i/35i2P6 



2/3 



(13) 



The thermal ICS gap is determined by those thermal-peak 
photons which has the maximum photon number density of 
the Planck spectrum at a certain temperature. It has a lower 
gap height hth 



hth = 2.7 X lO^P^/'^Pi 



-3/5^1/5j,-l/5 



(14) 



but lead to a larger Lorentz factor (72, th) for the secondary 
particles 



72, th 



3.7 X 103p-2/5p3/5^4/5j.V5 (15) 

(' Zhang etalJII997H) . 

These two CR and ICS gaps, which have relatively lower 
gap heights, would dominate the inner gap breakdown 
(I Zhang et all 1 1997bt iDu et al.ll2009l) . The pairs can also be 
abundantly generated by the primary particles escaped from 
the inner gap within a few neutron star radii, and they could 
have two major energy distributions due to the different types 
of primary particles. We will therefore use two different pair 
energy distributions for CR and ICS pairs to calculate the 
phase-averaged spectrum and phase-resolved spectra for the 
Crab pulsar. 

3. MODELING THE MULTI-WAVELENGTH LIGHT 
CURVES AND SPECTRA FOR THE CRAB PULSAR 

To explain the multi-wavelength light curves with nearly 
aligned peak phase for the Crab pulsar, we should obtain the 
high signal-to-noi se data, which are adopted from Figure 2 of 
dAbdo et alj|2010l) . We also reprocessed the Fermi 7-ray data 
to obtain the three 7-ray band light curves (see left panels of 
Figure O in the following steps: 

1. Limited by the timing solution for the Crab pulsafl 
from the Fermi Science Support Center (FSSC), we re- 
processed the original data observed from 2008 August 
4 to 2009 April 8. 

' http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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2. We selected photons of 0.1-300GeV in the "Diffuse" 
event class, within a radius of 2° of the Crab pulsar 
position (RA= 83.63°, DEC= 22.01°) and the zenith 
angle smaller than 105°. 

3. As done by lAbdo et alJ (120 lOl) . we used "fselect" to 
select photons of energy i?GeV within an angle of < 
max [6.68 — 1.76 log]^Q(_BGcv), 1.3] degrees from the 
pulsar position. 

4. We then obtained t he rotational phase for each photon 
using the tempo2 (iHobbs et al.l 120061) with the Fermi 
plug-in. 

5. Finally we obtained the multi-wavelength 7-ray light 
curves with 256 bins, as presented in Figure [3] (left pan- 
els). Two sharp phase-aligned peaks have a phase sep- 
aration of 5(f> ~ 0.4. 

An acceptable model should have reasonable input parame- 
ters (e.g., magnetic inclination angle a and viewing angle Q, 
and consistently produce multi-wavelength light curves with 
phase-aligned peaks and bridge emission and phase-resolved 
spectra for the Crab pulsar . 

3.1. Light Curves Modeling 

We adopted the method (see details in §3.1 of iDu et all 
I2OI IL including basic formulae and gap physics) to model 
the light curves. The key idea of modeling light curves is to 
project the radiation intensity of every spot on each open filed 
line (in either annular gap or core gap) to the "non-rotating" 
sky, with considerations of physical effects. Here the emissiv- 
ities are numerically assumed to facilitate calculations, they 
are however consisten t with the physic ally calculated spectra, 
as noted in Fig ure 8 of lDu et al.lfcOl ll) . Some model parame- 
ters for both the annular gap and core gap should be adjusted 
for the emission regions where the corresponding waveband 
emission are generated. The framework of the annular gap 
model as w ell as the details of coordinate system have been 
presented in iDu et akl (120101) . which can be used for simula- 
tion of the multi-wavelength hght curves of the Crab pulsar. 
We adopted the inclination angle of a = 45° and the viewing 
angle ( = 6 3° which were obtain ed from the Chandra X-ray 
torus fitting (iNg & Romanil2008l) . The modeUng methods are 
briefly delineated as follows. 

1. We first use the critical field line to separate the po- 
lar cap region into the annular and core gap regions. 
Then we use the so-called "open volume coordinates" 
('^ovcV's) to label the open field lines of the annular 
gap and core gap, respectively. Here rove is the nor- 
malized magnetic colatitude and tps is the magnetic az- 
imuthal. For the annular gap, we define the inner rim 
^ovc, AG = for the critical field lines and the outer 
rim rove. AG = 1 for the last open field lines; while for 
the core gap. We define the outer rim rove, CG = 1 for 
the critical field lines and the inner rim rove, CG = 
for the magnetic axis. We also divide both the an- 
nular gap (0 < rove, AG ^ 1) and the core gap 
(0.1 < rove, ee ^ 1) into 40 rings for calculation. 

2. When calculating the emissivities for modeling light 
curves, we postulate that the emissivities along one 
open field line have a Gaussian distribution rather than 



the frequently used assumption of the unif orm emis- 
sivity along an open field line dPvks & Ru dak 2003]; 
iFang & Zhang||2010l) for both the annular gap and the 
core gap. To justify this key assumption, we already 
used the realistic acceleration field to plot the flux 
against emission altitude al ong an op e n field line for the 
Vela, as shown in Fig. 8 of iDu et aP (1201 Ih . From that 
figure, one can clearly see that the flux is likely to have 
a Gaussian distribution against altitude near the peak 
position. The peak emissivities /p(6'p,V's) may fol- 
low another Gaussi an distribution against for a bunch 
of open field Unes (ICheng et al.1 120001: IPyks & RudakI 
I2OO3I: IFang & Zhangl2oTor As seen above, we use two 
different Gaussian distributions to describe the emissiv- 
ities on open field lines for both the annular gap and the 
core gap. The model parameters are adjusted to maxi- 
mally fit the observed multi-wavelength light curves. 

3. To derive the "photon sky-map" in the observer frame, 
we first calculate the emission direction of each emis- 
sion spot ne in the magnetic frame; then use a trans- 
formation matrix Ta to transform ne into ngpin in the 
spin frame; finally use an aberration matrix to trans- 
form iispin to nobscrvor = {hx, %, nj in the ob- 
server frame. Here 0o = arctan(ny/nx) and ( = 
arccos(nz/ -^Hx^ + riy^ + n^^) are the rotation phase 
(without retardation effect) of the emission spot with 
respect to the pulsar rotation axis and the viewing an- 
gle for a distant, nonrotating observer. The detailed 
calculations for the aberration effect can be found in 
iLee etalldlOlOh . 

4. We also add the phase shift 5(/)ret caused by the retarda- 
tion effect to the emission phase, that is = 0o — (50rot- 
There is no minus sign for (/jq because of our coordinate 
system different from other models. 

5. The "photon sky-map", defined by the binned projected 
emission intensities on the (0, Q plane, can be plotted 
in 256 bins (see middle panel of Figure [3]l. The corre- 
sponding light curves cut by a line of sight with a view- 
ing angle C — 63° are finally obtained. For the view- 
ing angle C = 63°, any magnetic inclination angles of 
a between 40° and 65° in the annular gap model can 
produce light curves with two peaks and a large peak 
separation similar to the observed ones. The emission 
from the single pole is favored for the Crab pulsar in 
our model. 

The modeled light curves from radio to TeV band are pre- 
sented in Figure [3] (black solid lines), and the key model pa- 
rameters are listed in Table [T] Emission from PI and P2 of 
multi-wavelength light curves originates from the annular gap 
region in the vicinity of the null charge surface, while bridge 
emission comes from the core gap region. We emphasize that 
all the multi-wavelength emission are originated from only 
one magnetic pole, our annular gap is therefore a single-pole 
magnetospheric model. 

Wit h well-coordinated efforts for pulsar radio timing pro- 
gram, lAbdo et all ( 120101) determined the phase lag between 
radio emission and 7-ray light curves. The first 7-ray peak 
comes earlier than the 1 .4 GHz radio pulse by a small phase 
of ~' 0.00853, but they are nearly phase-aligned. From light 
curve modehng, we can obtain the emission locations for each 
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Fig. 3. — Multi-wavelength (radio, X-ray and 7-ray) light curves for the Crab pulsar The observations are shown in the left p anels, and the data of RXTE, 
INTEGRAL, Comptel and Nancay ar e taken from Figure 2 of Abdo et al. (2010). The 25— 100 GeV light curve data are taken from lAleksic et alH 201 1), and the 
TeV (> 120 GeV) data are taken from lVERITAS Co llaboration et al. ( 201 1). The photon sky-map (middle panels) for an inclination angle and the coiTesponding 
modeled light curves (right panels) for a viewing angle of f = 63° are also presented, using the single-pole annular gap model. Our annular gap model can well 
explain the multi-wavelength light curves with phase-aligned peaks. 
(A color version of this figure is available in the online journal.) 
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Pulsar phase 

Fig. 4. — 7-ray and radio emission altitudes for the Crab pulsar with a 
viewing angle <^ = 63° in the annular gap model. The 7-ray peak and radio 
(1.4 GHz) emission altitudes ai'e nearly overlapped in the annular gap region, 
which leads to the phase-aligned peaks of the two energy band. While bridge 
emission are generated from the core gap region. 
(A color version of this figure is available in the online journal.) 



TABLE 1 

Model parameters for multi- wavelength light curves of the 
Crab pulsar 



Band 


K 


A 


€ 




<ye. A 


CTC 


"■e, ci 


o"e, C2 


25-100 GeV 


0.65 


0.6 


0.8 


0.4 


0.0045 


0.25 


0.0023 


0.0035 


>I20GeV 


0.68 


0.85 


0.8 


0.3 


0.0035 


0.25 


0.0025 


0.0032 


2-16 keV 


0.50 


0.8 


1. 17 


0.5 


0.005 


0.3 


0.0022 


0.0035 


100-200 keV 


0.50 


0.8 


0.8 


0.6 


0.006 


0.35 


0.0035 


0.0045 


0.75-30 MeV 


0.52 


0.8 


0.7 


0.5 


0.004 


0.3 


0.005 


0.006 


0.1-0.3 GeV 


0.53 


0.8 


0.6 


0.5 


0.004 


0.15 


0.0012 


0.0012 


0.3-1.0 GeV 


0.54 


0.8 


0.7 


0.5 


0.005 


0.15 


0.0014 


0.0014 


> 1.0 GeV 


0.57 


0.8 


0.5 


0.5 


0.005 


0.35 


0.0035 


0.0035 



Note. K and A are two geometry parameters to determine the peak altitude in 
the annular gap; e is a parameter for the peak altitude in the core gap; (ja 
and (7(3 are length scales for the emission region on each open field fine in 
the annular gap and the core gap in units of i?LC. respectively; ag a is the 
transverse bunch scale for field fines in the annular gap; 0-9 ci and erg ^2 
are the bunch scale for field lines of —180° < i/Jb < 90° and 
90° < ipB < 180° in the core gap, respectivel y. The detailed description of 
these symbols can be found in IDu et alJ20rih . 



band. The result of 7-ray and radio emission altitudes for the 
Crab pulsar is shown in Figure H) The region for the radio 
emission is mainly located at a altitude of ~ 0.5i?LC on some 
certain filed lines, for PI, the magnetic azimuthal ip is in the 
range of -98° to -96°, while 106° to 109° for P2. It is gen- 
erated in an intermediate-altitude annular gap region, which 
might be due to the coherence condition and propagation ef- 
fects. It is found that the positions of both 7-ray (X-ray) peak 
and radio peak (1.4 GHz) are overlapped, this leads to the 
nearly phase-aligned pulse peaks of multi-wavelength emis- 
sion, except for several GHz radio emission due possibly to 
plasma propagation effects. Nevertheless, the 7-ray emission 
altitudes are above th e lower bound of the height determined 
by 7 — i? absorption (iLee et al.ll2010l) . Based on our model, 
not all 7-ray pulsars can be detected in the radio band, and not 
all radio pulsars can have a 7-ray beam with sufficiently high 
flux towards us. The beam shapes and intensities of 7— ray 



and radio can evolve with pulsar ages. 

3.2. Multi-wavelength Spectra for the Crab Pulsar 

In this section, we will use the annular gap model to cal- 
culate the multi-wavelength phase-averaged s pectrum and 
pha se-resolved spectra for the Crab pulsar iKuiper et al.l 
(I2OOL) achieved seven band phase-resolved spectra (LWl, PI, 
TWl, Bridge, LW2, P2, TW2) and the phase-averaged spec- 
trum of the Crab pulsar with the EGRET 7-ray data. We will 
also add the new high-quality Fermi, MAGIC and VERITAS 
data to the total phase-averaged spectrum. The 7-ray emis- 
sion is believed to be originated from t h e curvature radiation 
of primary partic les (iTang et al.l 120081: iHarding et all 120081: 
iMeng et al1l2008l) . which generally gives a super-exponential 
power-law spectrum with the cut-off energy of a few GeV. 
However, the TeV (20 to 400 GeV) emission mechanism re- 
mains a mystery, and it requires a global phase-averaged spec- 
trum fitting to resolve this problem. 

We first discuss the particle dynamics in the annular gap. 
After exploring the formation mechanism of acceleration 
electric field, the relevant dynamic parameters: acceleration 
electric field E\\, Lorentz factor of primary particles 7p, char- 
acteristic energy of curvature radiation emitted from primary 
particles Ec and escape photon energy Ecs are calculated for 
two field lines where PI and P2 are mainly originated. The 
calculated results are shown in Figure |5] As introduced in 
§ 2.2, the derived acceleration electric field E\\ is quite huge 
in the inner region of annular gap below the altitude of i?ii,n, 
which leads to the generation of numerous pairs via 7— B ab- 
sorption effect. The Lorentz factor of primary particles 7p is 
derived from the curvature radiation reaction using the Equa- 
tion (O. However, the actual Lorentz factor of the primary 
particles is smaller than the one shown in Figure |5] if other 
energy loss mechanisms (e.g. ICS loss rate) as well as rela- 
tivistic and pair screening effects are taken into account. Ec 
denotes for the characteristic energy of curvature radiation 
emitted from primary particles, derived by the Equations Q 
and ( [Tol l. Ecs, due to the 7— B absorption effect based on the 
Equation ( |26] l. is the escape (maximum) photon energy for a 
certain altitude. 

We will further consider the emission mechanisms for the 
Crab pulsar. The photon number spectrum of the synchrotron 
radiation is given by 



d^i. syn 



V3e^B{r) sin ip{r) 1 

TTleC^h E^ 



G{x), (16) 



where 95 (r) is the pitch angle at a distance of r on an open 
field line, h is the Planck constant, G{x) = x Ks {£,)d^, 

K5 is the modified Bessel function with an order of 5 /3, x = 

3 

E^/Esyn,c{r) and 



= 4.3 X 10^ hB{r)j'^ sm(p{r) 



(17) 



is the critical synchrotron photon energy. The pitch angle 
(p{r) of relativistic primary particles flowing along a mag- 
netic field line could be small, but it cannot be neglected for 
synchrotron radiation. While the pitch angle (p{r) of pairs 
increases due to the cyclotron resonant absorption of the low- 
energy photons (Harding et al. 2008) and it varies with the 
emission altitudes. The mean pitch angles of the two types 
of pairs are different owing to the effect of cyclotron resonant 
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Fig. 5. — Four dynamic parameters along two open field lines where PI (top panel) and P2 (bottom panel) are mainly originated. The acceleration electric field 
(solid fine) is quite huge in the inner region of annular gap below the altitude of Rum, which results in the generation of numerous pairs. The Lorentz factor 
of primary particles 7p (dotted line) is derived from the curvature radiation reaction. Eq (dashed line) denotes for the characteristic energy of curvature radiation 
emitted from primary particles. Ecu (dot-dash line) is the escape (maximum) photon energy for a certain altitude due to the 7— B absorption. 



absorption for different particles with different Lorentz fac- 
tors. The synchrotron radiation from pairs play an important 
role in X-ray band to soft 7-ray band (e.g. < 0.02 GeV). 

The photon number spectrum of the curvature radiation can 
be given by 



\dE^dt 



6 7 



hp(r)Ej 



(18) 



where p(r) is the curvature radius at r and Equation ( fTOl i 
s hows the critical curva t ure ph oton energy. 

iBlumenthal & Gouldl ( 119701) presented an analytic formula 
for the photon spectrum of the inverse Compton scattered 
photons per elect ron in the case of extreme Nishia- Klein limit 
and then lTang et a l. (2008) gave a simplified form, i.e.. 



3o-TC nsyn(£, r) +nx(£,r) 
472 e 



2q\nq+{l + 2q){l-q) + 



(rg)^(i-g) 

2(1 + rg) 



de, (19) 



where q = Ei/r{l-Ei), T ^ A-ie/raeC^ and £^1 = E^/E^. 
£ is the energy of soft photons for scattering, and ei and £2 
are the minimum and maximum energy of the soft photons 
for integration, respectively. We choose the values of £1 and 
£2 to fulfill the condition of Ei < 1. The lower limit £1 in the 
Equation(fT9ll is chosen to be around 1 eV, and the upper limit 
£2 is adjusted artificially to make a quick convergence of the 
Equation(fT9ll. 

There are two possible sources of soft photons for the in- 
verse Compton scattering, one is the thermal photons, and the 
other is synchrotron photons. The thermal photons, generated 
by the stellar surface with a typical surface temperature T, 
is an important source for Compton scattering of the primary 
particles and secondary particles. The number density of soft 



photons is given by 



nx(£,r) 



R 

■K'^{hcf exp{e/kT) - 1 V~ 



(20) 



where k is the Boltzmann constant. The temperature T is 
taken to be 2 x 10^ K in our calculations. 

The other source of soft photons arises from synchrotron ra- 
diation of pairs. Owing to the quite abundant soft synchrotron 
photons, the scattering of this kind of photons is more signif- 
icant at higher altitudes near the light cylinder Although the 
synchrotron radiation spectrum from a single particles with a 
Lorentz factor of 7 is maintained to a wide energy band (for 
example, soft X-ray band to 7-band), it is more likely to be 
a spectral line in fact, which is concentrated on the critical 
energy. Thus we can rewrite the synchrotron spectral power 
using the total energy loss rate, i.e., 

P^^syn ^ EsynS{s ~ Esyn. c) ^ ^(TtcJ7b7^^(£ - Ec) , (21) 

where ctt is the Thompson scattering cross section, Ub is the 
energy density of the local magnetic field, and (5(£ — £c) is a 
Delta function. 

Supposing pairs follow a power-law spectrum neij) = 
no7~'' with a particle energy spectral index of s, we can ob- 
tain the synchrotron emissivity in a simple form 

(22) 

where z^l is the Larmor frequency of an electron in the mag- 
netic field, no is a constant and Ail(r) is the solid angle of 
the beam of synchrotron photons and is estimated as 



27r /"fiir) 



JO 



)sm 



TTip^ir). (23) 
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TABLE 2 

Best fit parameters for the calculated phase-averaged spectrum of the Crab pulsar 



c 


■yCR 
! mm, pair 


-7CR 

'max, pair 


pair 


-vies 

'mm, pair 


-vies 

'max, pair 


"pair 


pri 
' min 


pri 

/max V^pri 


45° 63° 


2.5 X 10^ 


2.18 X 10'' 


0.0092 


6.0 X 10^ 


1.19 X 10^ 


0.0074 


1.0 X lO'* 


2.78 X lO'^ 0.000098 



TABLE 3 

Best fit parameters for the calculated phase-resolved spectra of the Crab pulsar 



Phase band 


'mm, pair 


-yea 

' max, pair 


«CR 

" pair 


-yies 

' mm, pair 


-yies 

'max, pair 


<pics 

" pair 


7P" 

'mm 


pri 
Tmax 






LWl 


2.5 X 10^ 


1.2 X 10'' 


0.009 








1.0 X 10" 


1.76 X 10 ' 


0.00005 


1.82 


PI 


2.5 X 10^ 


0.9 X lO" 


0.0078 


4.6 X 10^ 


0.86 X 10^ 


0.0071 


1.0 X 10" 


2.47 X lO'^ 


0.00009 


3.03 


TWl 


2.5 X 10^ 


0.95 X 10* 


0.009 


3.0 X 10^ 


0.68 X 10^ 


0.0079 


1.0 X 10^ 


2.78 X lO'^ 


0.00008 


3.16 


Bridge 


3.0 X 10^ 


0.71 X 10* 


0.0085 


5.0 X 10^ 


0.55 X 10^ 


0.0070 


1.0 X 10" 


2.55 X lO'^ 


0.00003 


2.53 


LW2 


3.5 X 10^ 


1.25 X 10" 


0.008 


5.0 X 10^ 


0.52 X 10^ 


0.0079 


1.0 X 10" 


2.50 X lO'^ 


0.00005 


1.78 


P2 


5.0 X 10^ 


1.22 X lO* 


0.0078 


5.0 X 10^ 


0.82 X 10^ 


0.0085 


1.0 X 10*^ 


2.97 X lO'^ 


0.000072 


3.81 


TW2 


5.0 X 10^ 


1.3 X 10" 


0.008 








1.0 X 10'^ 


1.82 X lO'^ 


0.00007 


2.78 
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FiG- 6. — The modeled phase-averaged spectrum for the Crab pulsar. The 
calculated total spectrum (thick black solid line) is obtained from the sum 
of synchrotron radiation and ICS from two kinds of pairs and curvature ra- 
diation and synchrotron radiation from primary particles. It is found that 
the curvature radiation and synchrotron radiation from primary pai'ticles is 
mainly contributed to 7-ray band (20 MeV to 20GeV); synchrotron radia- 
tion from CR-induced pairs and ICS-induced pairs dominates the X-ray band 
and soft 7-ray band (lOOeV to lOMeV). ICS from the pairs contributes to 
hard TeV 7 -rav band 20 GeV to 400 GeV). The data (solid circles) are 
taken from IKuiper et alj j200l|). The Fermi 7-ray data is plotted in pur- 
ple circles. The 7-ray spectral data of VERITAS (square) are taken from 
[VERITAS Collaboration et al. ( 20 1 1 ) . and the MAGIC data (big hollow cir- 
cle) are taken from MAGIC Collaboration et al.l 120121) . 
(A color version of this figure is available in the online journal.) 

Therefore we use an approximate formula to facilitate ICS 
calculations when we consider the soft seed synchrotron pho- 
tons. The number density of the soft synchrotron photons, 



syn(e, r) can be given by 
Aflr . 



2(tt?7b?t-o?' 
3e£ 



(24) 



where = hi^ is the Lar mor energ y. 

Following the method of lDu et aiT(l201 lb . we divide the an- 
nular gap region into 40 rings and 360 equal intervals in the 
magnetic azimuth, i.e. in total 40x360 small magnetic tubes. 



A small magnetic tube has an area ^0 on the neutron star 
surface. From Equation (|2]), the number density of primary 
particles at a altitude r is n(r) = cos^out- The cross- 

section area of the magnetic tube at r is A{r) = BqAq/ B{r). 
The flowing particle number at r in the magnetic tube is 



AN (r) ^ Air) As— ^ cos Cant, 

ZTTCe 



(25) 



here As is the arc length along the field. The energy spectrum 
^pri = dN/dj of the accelerated primary particles is not well 
understood in the first physical principle. Here we assume the 
primary particles in the magnetic tube to follow a power-law 
energy distribution diVTd-y = No-f^ with an index of F = 
—2.2. Using Equation ( |25] |, iVo can be derived by integration 
the equatio n above. 

Harding et al.l (120081) have assumed a broken power-law 
distribution for pairs with indexes of —2.0 and —2.8 [see their 
Equation (47)]. But in our model, we postulate that the two 
types of relativistic pairs follow two different power-law en- 
ergy distribution as noted in § 2.3, i.e. 



CR 
max' 



^-L / ' 'mm ^ ' ^ '- 

-.2 yes ^ ICS 

■ ^ / ' 'mm ^ / ^ / max' 



where si = 2.45 and S2 = 2.6 are the spectral index, Ci and 
C2 are two coefficients, 7^^ , 7,?^^, and j}^^^ are lower 
limits and upper limits of the Lorentz factors for the CR and 
ICS pairs. 

The pairs can be generated with a large multiplicity (lO'^ — 
10'') via the 7 — S process in the lower regions of the an- 
nular gap and the core gap near the neutron star surface. The 
pitch angle of pairs increases due to the cyclotron resonant ab- 
sorption of the low-energy photons. The mean pitch angle of 
secondary particles is considerable owing to the effect of cy- 
clotron resonant absorption. The synchrotron radiation from 
secondary particles has some contributions to the low-energy 
7-ray emission, e.g. < 0.02 GeV. 

We also had an analytic al formula of op tical depth t^^b 
due to the 7— B absorption (iLee et al.ll2OI0h 



(r) 



1.55 X lOV 



2.76 X io6r5/2pi/2 



)' (26) 



^7 J->{>, LAJ^l- J^J 

here is in units of MeV, Bq, 12 is the surface magnetic field 
in units of 10^^ G. We found that the photons of the Crab pul- 
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Fig. 7. — Similar as Figure |6] but for modeled seven phase-band phase-resolved spectra of the Crab pulsar. TWl a nd LW2 are plotted twice. The MAGIC 
(50—400 GeV) spectral data are available only for PI and P2, which are taken from lMAGIC Collaboration et alj 42012I) . whereas the Fermi and VERITAS data 
are not included here. 

(A color version of this figure is available in the online journal.) 



sar with an energy of < 100 GeV always satisfy the condition 
of T-y-B ^ 1 if the emission altitude is greater than a few 
hundred kilometers. Thus the final multi-wavelength spec- 
trum emitted by the primary particles and secondary particles 
can be calculated by 



\dE^dt 

r^-B (r) 



dE-ydt , 

7 / syn 



(27) 



where A17off is the effective solid angle of the emission beam, 
and £> = 2 kpc is the distance from the Crab pulsar to the 
Earth. 

There are actually six main spectral components for the pul- 



sar total spectrum. Based on our calculations, the synchrotron 
and curvature radiation from primary particles and the syn- 
chrotron radiation and ICS from secondary particles are re- 
quired to calculate the phase-averaged and phase-resolved 
spectra, while ICS from primary particles and CR from pairs 
can be ignored for the Crab pulsar. When we calculate the 
SR spectra, we carefully deal with the modified Bessel func- 
tion with for unprecedent accuracies. While calculating ICS 
spectrum, as stated above, we use Equations ( fT9] l and ( l20b and 
( l24l i to reduce the computation time. 

To accelerate our computations, we use t he method of " av- 
eraged emission-altitude" (as introduced in lDu et al]l201 II) to 
calculate the four contributive spectral components for the 
Crab pulsar, i.e., synchrotron and curvature radiation from 
primary particles, and ICS and synchrotron radiation from 
pairs. We first obtain the emission altitudes from Figure 
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(|4|i for each phase band. For instance, the emission altitude 
is about 0.54 i?LC on the field line of a magnetic azimuth 
ip = —101° for PI; while for P2, the emission altitude is 
0.49 i?LC on the field line of V = 109°. Then we com- 
pute their acceleration electric field and potential ^P, and 
adjust the minimum and maximum Lorentz factor for both 
primary particles and pairs (7j'J"„ and 7,p"x Imit ^^'^ 
7max)' the pitch angle e(r) and the 7-ray beam angle Ailcff to 
fit the multi-wavelength phase-averaged spectrum and phase- 
resolved spectra for the Crab pulsar 

We fitted the phase-averaged spectrum and seven phase- 
band phase-resolved spectra of the Crab pulsar, and the 
results are shown in Figure |6] and Figure |7] respectively. 
For phase-averaged s pectrum, we basica lly used the multi- 
wavelength data from lKuiper et al.l (120011) and combined with 
latest Fermi, MAGIC and VERITAS 7-ray spectral data, 
whereas for phase-re solved spectra, we o nly used the multi- 
wavelength data from lKuiper et al.l ( 1200 Ih . The best fit param- 
eters for phase-averaged spectrum and phase-resolved spectra 
are listed in Table |2] and Table [5] respectively. From spec- 
tra fitting, we found that the calculated 7-ray spectra are not 
sensitive to 7^^"^, but quite sensitive to j^l^ which is cho- 
sen below the value obtained from the balance of curvature 
radiation and radiation reaction shown in Figure |5] The solid 
angle Af^cff was always assumed to be 1 by many authors for 
simplicity. We adjusted it as a free parameter around 1 for 
different phases. 

We found that multi-wavelength emission from the phase 
bands of PI, P2 and bridge contribute significantly to the to- 
tal phase-averaged spectrum. The phase-averaged spectrum 
and phase-resolved spectra are decomposed into four spectral 
components. Curvature radiation and synchrotron radiation 
from primary particles is the main origin of the observed 7- 
ray emission (lOMeV to ^ 20GeV), synchrotron radiation 
from CR-induced pairs and ICS-induced pairs dominates the 
X-ray band and soft 7-ray band. ICS from the pairs con- 
tributes significantly up to TeV band, while ICS from both 
primary particles and curvature radiation from pairs can be 
neglected for the spectrum fitting. Owing to the larger emis- 
sion altitudes (which leads to lower acceleration electric field) 
for LWl and TW2, the lorentz factors 7 are very low which 
results from lower acceleration electric field, curvature radia- 
tion and synchrotron radiation have little contributions to the 
LWl and TW2 phase band spectra. The TeV emission of ICS 
from pairs can be also found for PI, bridge and P2 in our cal- 
culated phase-resolved spectra, which are consistent with the 
modeled TeV light curve. Two types of pairs, CR-induced and 
ICS -induced, could be therefore confirmed by the spectra of 
the Crab pulsar, these may be also the origin of the two types 
of wind pairs in Crab Nebula. 

4. CONCLUSIONS AND DISCUSSIONS 

Owing to its strong multi-wavelength emission, the famous 
Crab pulsar is a crucial astrophysical object to distinguish the 
emission mechanisms from different magnetospheric models. 
In this paper, we calculated radio. X-ray, 7-ray and TeV light 
curves, phase-averaged spectrum and phase-resolved spectra 
in the framework of the 3D annular gap and core gap model 
with reasonable emission-geometry parameters {a = 45° and 
C = 63°). It is found that the electric field in the annular gap 
is huge (~ 10^*^ eV) in the several tens of neutron star radii 
and vanishes near the region of i?LC- The primary particles 
are accelerated to ultra-relativistic energies, and produce nu- 



merous secondary particles (CR and ICS pairs) in the inner 
region of the annular gap via 7— B process. The pulsed emis- 
sion of radio. X-ray and 7-ray are generated from the emission 
of primary particles or secondary particles (pairs) with differ- 
ent emission mechanisms in the nearly similar region of the 
annular gap (or core gap) in only one pole, this leads to the 
"phase-aligned" multi-wavelength light curves. The emission 
of PI and P2 is originated from the annular gap region near 
the null charge surface, while the emission of bridge is mainly 
originated from the core gap region. 

Assuming that power-law energy distributions of primary 
particles and two types of pairs, the phase-averaged spec- 
trum and phase-resolved spectra of the Crab pulsar are well 
produced by mainly four components: synchrotron radiation 
from CR-induced and ICS-induced pairs dominates the X- 
ray band to soft 7-ray band (lOOeV to lOMeV); curvature 
radiation and synchrotron radiation from the primary parti- 
cles mainly contribute to 7-ray band (lOMeV to ~ 20GeV); 
ICS from the pairs significantly contributes to the TeV 7-ray 
band (100 GeV to 400 GeV). Our fitted phase-averaged spec- 
trum and phase-resolved spectra have similar tendency vary- 
ing with the phot on energy and are basically consistent with 
outer gap model ( | Tang et aP l2008h and the slot gap model 
(iHarding et al.ll2008l) ' at soft X-ray to a few tens of GeV band, 
but quite different in > 20 GeV band. This is mainly due 
to the additional spectral component of ICS from pairs. From 
the multi-wavelength spectral fitting we emphasize that curva- 
ture radiation alone emitted from the primary particles cannot 
explain the TeV band (25 to 400 GeV) emission of the phase- 
averaged spectrum, ICS from pairs significantly contribute to 
this > 20 GeV band. In addition, two types of pairs are gen- 
erated in the magnetosphere, and they may be also the origin 
of the two types of wind pairs in Crab Nebula ( Abdo et aT] 
I20T0I) . 

Radio emission ( 1 .4 GHz) of the Crab pulsar is originated 
from a narrow and high-altitude region with a similar location 
of the 7-ray emission, which leads to the phase-aligned peaks. 
Our model for radio photon sky-map is patch-like, however 
the detailed emission mechanism for radio emission is needed 
to further studied. 
The popular outer gap ( exc ept the versions of 



Romani & Yadi garoglul 1119951): |Cheng et al.] (120001) : 
Romani & Watte rsI ( 201 Oh) and slot ga p model s are two- 



pole models (Tang et al.''2()08l :IZhang & Liii2009i:lLi & ZhangI 
[2OIO; Harding et al. 2008^ To model the observed light 
curves and spectra for the Crab pulsar, they require the 
emission from the both magnetic poles, which result from 
larger magnetic inclination angle a or larger viewing angle 
C,. However, our annular gap model is an intermediate 
emission-altitude and single-pole mo del with reasonable a 
and ( from the X-ray torus fitting ( iNg & Romanil |2008l) . 
Unfortunately, the important paramete rs a and C for pulsar 
emission geometry are uncertain so far iNg & Romanil ( 120081) 
can only give a reliable viewing angle ( for some young 
pulsars which have X-ray torus configurations, and then 
combine with the radio rotating vector model (RVM) to 
obtain the inclination angle a using the radio polarization 
angle (PA) fitting. The simple RVM model is only based on 
the geometry at a certain low altitude for an assumed circular 
emission beam, and the propagation effects that can change 
the polarization states that had been already ignored by the 
RVM model. The derived a by this method is therefore 
debatable. A better method is strongly desired to obtain the 
convincing values of a and C- 
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To well explain the multi-wavelength pulsed emission 
from pulsars, the detailed magnetic field configuration and 
3D global acceleration electric field distribution with proper 
boundary conditions for the annular gap and the core gap 
should be carefully studied. Unfortunately , these two physical 
aspect s are not fully understood. Recnetlv. lRomani & Walters! 
(120101) studied pulsar light curves with magnetosphere beam- 
ing models and found that outer gap model and approximat- 
ing force-free dipole field were preferred a t their high statisti- 
cal significance. However. iHarding et al.l (1201 ll) also studied 
high-energy pulse profiles (e.g. the Vela pulsar) using both 
retarded vacuum dipole and force-free field geometry. They 
found that slot gap model with vacuum dipole was more fa- 
vorable. Therefore, the subject of pulsar magnetic field con- 
figuration is still debatable. In addition, the problem of 3D 
acceleration field with the general relativistic effect and pair 
screening effect is more complicated, although many efforts 
have been paid. We just derived the ID (actually 2D) contin- 
uous solution for the acceleration Possion equation, and the 
general relativistic and pair screening effects have not been 
taken into account in our annular gap model at present. This 
is our first step to establish our model picture, and will benefit 



further 3D complicated physical studies with considerations 
of related effects. We emphasize that some simplified hypoth- 
esis considering qualitative physical effects have been used 
in our model to study pulsar light curves and spectra. This 
can give us insightful enlightenments to improve our knowl- 
edge of pulsar radiation physics. We will further improve our 
model to give more precise modeled light curves, especially 
for the phases of LWl and TW2. 

In sum, the multi-wavelength emission from the Crab pul- 
sar can be well explained in the annular gap and core gap 
model, and this is also done for the Vela pulsar (iDu et al.l 
lIoTlh . Our model is a promising model to unveil the multi- 
wavelength pulsed emission from 7-ray pulsars. 



The authors are very grateful to the referee for the insightful 
and constructive comments. We thank both the pulsar groups 
of NAOC and of Peking University for useful conversations. 
The authors are supported by NSFC (10821061, 10573002, 
10778611, 10773016, 11073030 and 10833003) and the Key 
Grant Project of Chinese Ministry of Education (305001). 

Facilities: Fermi (LAT), MAGIC, VERITAS 



REFERENCES 



Abdo, A. A., et al. 2010, ApJ, 708, 1254 

Aleksic, J., Alvarez, E. A., Antonelli, L. A., et al. 201 1, ApJ, 742, 43 
Blumenthal. G. R., & Gould, R. J. 1970, Reviews of Modem Physics, 42, 
237 

Cheng, K. S., Ho, C, & Ruderman, M. 1986a, ApJ, 300, 500 
Cheng, K. S., Ho, C, & Ruderman, M. 1986b, ApJ, 300, 522 
Cheng, K. S., Ruderman, M., & Zhang, L. 2000, ApJ, 537, 964 
Daugherty, J. K., & Harding, A. K. 1994, ApJ, 429, 325 
Daugherty, J. K., & Harding, A. K. 1996, ApJ, 458, 278 
Du, Y. J., Xu, R. X., Qiao, G. J., & Han, J. L. 2009, MNRAS, 399, 1587 
Du, Y. J., Qiao, G. J., Han, J. L., Lee, K. J., Xu, R. X. 2010, MNRAS, 406, 
2671 

Du, Y. J., Han, J. L., Qiao, G. J., & Chou, C. K. 201 1, ApJ, 731, 2 
Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201 
Erber, T. 1966, Reviews of Modern Physics, 38, 626 
Fang, J., & Zhang, L. 2010, ApJ, 709, 605 
Goldreich, P, & Julian, W. H. 1969, ApJ, 157, 869 
Harding, A. K., Stem, J. V., Dyks, J., & Frackowiak, M. 2008, ApJ, 680, 
1378 

Harding, A. K., DeCesar, M. E., Miller, M. C, Kalapotharakos, C, & 

Contopoulos, I. 201 1, arXivil 1 1 1.0828 
Hirotani, K. 2008, ApJ, 688, L25 

Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, MNRAS, 369, 655 
Kuiper, L., Hermsen, W., Cusumano, G., Diehl, R., Schonfelder, V., Strong, 

A., Bennett, K., & McConnell, M. L. 2001, A&A, 378, 918 
Lee, K. J., Du, Y. J., Wang, H. G., Qiao, G. J., Xu, R. X., & Han, J. L. 2010, 

MNRAS, 405, 2103 
Li, X., & Zhang, L. 2010, ApJ, 725, 2225 
Lin, G. F, & Zhang, L. 2009, ApJ, 699, 1711 



MAGIC Collaboration, Aleksic, J., Alvarez, E. A., et al. 2011, 
farXiv: 1109.6124 

Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 
1993 

Meng, Y, Zhang, L., & Jiang, Z. J. 2008, ApJ, 688, 1250 
Muslimov, A. G., & Harding, A. K. 2003, ApJ, 588, 430 
Muslimov, A. G., & Harding, A. K. 2004, ApJ, 606, 1 143 
Ng, C.-Y, & Romani, R. W. 2008, ApJ, 673, 41 1 
Qiao, G. J., & Lin, W. P 1998, A&A, 333, 172 

Qiao, G. J., Lee, K. J., Wang, H. G., Xu, R. X., & Han, J. L. 2004a, ApJL, 
606, L49 

Qiao, G. J., Lee, K. J., Zhang, B., Xu, R. X., & Wang, H. G. 2004b, ApJL, 
616, L127 

Qiao, G. J., Lee, K. J., Zhang, B., Wang, H. G., & Xu, R. X. 2007, Chinese 

Joumal of Astronomy and Astrophysics, 7, 496 
Romani, R. W., & Yadigaroglu, I.-A. 1995, ApJ, 438, 314 
Romani, R. W., & Watters, K. P 2010, ApJ, 714, 810 
Ruderman, M. A., & Sutherland, P G. 1975, ApJ, 196, 51 
Tang, A. P S., Takata, J., Jia, J. J., & Cheng, K. S. 2008, ApJ, 676, 562 
VERITAS Collaboration, Aliu, E., Arlen, T, et al. 201 1, Science, 334, 69 
Zhang, B., Qiao, G. J., Lin, W. P, & Han, J. L. 1997a, ApJ, 478, 313 
Zhang, B., Qiao, G. J., & Han, J. L. 1997b, ApJ, 491, 891 
Zhang, L., & Cheng, K. S. 1997, ApJ, 487, 370 

Zhang, L., Cheng, K. S., Jiang, Z. J., & Leung, P 2004, ApJ, 604, 317 
Zhang, L., Fang, J., & Chen, S. B. 2007, ApJ, 666, 1165 

Zhang, L., & Li, X. 2009, ApJ, 707, L169 



